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Craps et al. recently proposed to study a simple cosmological background with a null 
singularity [0, this cosmology admits a matrix model description, thus lends itself to a 
rigorous study. An earlier example of cosmology with a null singularity is proposed in 
and is later studied by many authors , they find that this singularity is highly unstable 
due to gravitational back-reaction. However, the model of Craps et al. seems to avoid this 
problem. 

The construction of the model in is simple. One starts in type IIA string theory 
with a fiat string metric and a linear dilaton background (j) = —Qx^, where is a light- 
like coordinate. The linear dilaton does not require modifying the critical dimension, since 
its linearity is along a null direction. Although the string metric is fiat, the Einstein metric 
is nontrivial: 

ds^ = e^^^/^{-2dx+dx- + {dx'f). (1) 

For a positive Q, the metric contracts to a singularity at x~^ = — oo, this is actually 
a curvature singularity. The corresponding 11 dimensional M theory metric also has a 
singularity at x~^ = oo. The authors of |jl[] show that the singularity at x~^ = — oo lies in a 
finite geodesic distance, while the other singularity is at infinite distance. This background 
preserves half of total 32 supersymmetries, one expects that there is a control over the null 
singularity. 

Since the string metric is fiat, this time-dependent background admits a matrix string 
description, in which the Yang-Mills coupling constant is time-dependent. In fact, the 
authors of Q show that this Yang-Mills theory can be regarded as one with a constant 
coupling on a world sheet with a time-dependent metric. 

It goes without saying that this is an important observation, and it may provide the 
first example in which we can study a time-dependent background in a controlled fashion. 
Thus, it is interesting to ask whether this model is unique, or it has many cousins. In this 
note we will see that indeed there is a large class of such models. 

We shall limit ourselves in M theory in this note. Consider metric 

ds^ = e^"^^ {-2dx+dx- + {dx')^) + e^^^^ [dx")^, (2) 

where i = 2, ... 10 — (i, a = 11 — (i, . . . 10, namely, there are d coordinates x"", the total 
dimensions of spacetime is 11. The metric of [|1| is given by taking d = 1, a = Q/3 and 
/3 = —2(5/3. X^^ is taken as the M theory circle, thus the metric reduces to the IIA fiat 
string metric with a linear dilaton background. 
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To find a solution, we use an orthonormal basis 

e± = e""^rfa;±, = e"""^ dx\ e" = e^'^rfx". (3) 
The non-vanishing components of spin connection are 

u+- = ae-""^e+, u,+ = ae-""^e\ cc;„+ = /3e-""^e". (4) 
The non-vanishing curvature 2-forms are 



iti_i_ = a e e A e , 



(5) 



= /3(2a - /3)e-2"^ e" A e+. 
The only non-vanishing component of the Ricci tensor is R++ given by 

R++ = [(9 - d)a^ + d(3{2a - /3)]e-2"^^ , (6) 

and the Einstein equation R++ = has two solutions 

/3 = (1±-^K (7) 

For (i = 1, choosing the minus sign in (|^) reproduces the background considered in 

We now show that the background given by eqs.(|l) and (|^) not only is a solution, but 
also preserves half of super symmetries. The only interesting SUSY transformation is that 
for the gravitino = -D^e. The components of spin connection of interest are uj+, uji 
and Ua- 

ojj^ — —2a — 2a7~7"'", Ui = 207*7"'" 
Ua = 2/3e(^-")^^7"7+. 

Since 



(8) 



D,] = - 2a)e(^-)^^7 V, (9) 

The compatibility condition for the constraints D^e = is 7"'"e = 0. This condition 
eliminates half of the components in e. The constraints Die = Dae = are solved if e is 
independent of and x". Finally, -D+e = is solved if e = exp(iax"'")77 for a constant rj. 
We conclude that the unbroken supersymmetry is parametrized by rj with the constraint 
7"'"r7 = 0. Thus, just like the original background of [Q, our more general metric also 
preserves 16 supersymmtries, this should be enough to guarantee that our proposed matrix 
model to be described shortly is a valid description of dynamics over this background. 
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In the case d = 9, (3 = 0, 2a, although [-D+, Da] = 0, but [-D+, Di] = —| 0^7*7+, again 
we obtain the same unbroken supersymmetries. 

In fact, d^) is only a special case of a more general class of solutions preserving half 
of supersymmetries (to be discussed in the end of this note). Before we study the matrix 
model for (^, let us discuss the geometric properties of (0). Apparently, the sign of a and 
its absolute value can be changed by changing and x~ , so we always assume a > 0. As 
x^ —>■ —00, the factor e^"^^ approaches zero, the transverse dimensions shrink to zero 
size, this is the big bang point with regard to these coordinates. This singularity locates 
a finite distance away in view of the affine parameter, indeed, let dX^ = e'^"^^ dx^ , the 
first term in the metric (|]) becomes 

-2dX+dx-, (10) 

thus is the affine parameter for the null geodesic x~ =const. The big bang singularity 
occurs at X^ = 0. 

The geometry of other transverse dimensions depends on the choice /3 in (^. We 
name the choice (3= {1 + 3/y/d) a case 1, and the choice /3 = (1 — 3/^/d)a case 2. For case 
1, /3 > 0, so all x"" shrink at x~^ = —00. For case 2, /9 < except for the d = 9 case, thus 
all x"" get to infinity at x~^ = —00 and shrink to zero size at x^ = 00. d = 9 is special, in 
this case, all 9 dimensions do not evolve in time, and one can redefine x^ such that the 
metric is Minkowski. Supersymmetry is also enhanced, since there is no longer constraint 
7+e = following from (^. 

The model studied in is a special example of case 2 when d = 1. As in we can 
compactify x^^ on a circle to obtain type IIA string theory. The string metric and the 
dilaton (p are related to the M theory metric through 

ds^^e-'^/'dsl + e^'^/\dx'^f, (11) 

we obtain 

(P = ^/3x+, ds^t = e(2"+/5)^^[-2c/x+dx- + {dx^f] + e^^""^ {dx^f . (12) 
The lOD Einstein metric reads 

dsl = e(2«+/3/4)x+ ^_2dx+dx- + [dx'f] + e(9^/4)-^ {dx'^f. (13) 
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Since the string coupling "constant" Qs = exp((^) = exp(|/3a;"'"), for case 1, strings are 
weakly coupled at big bang, while all transverse dimensions are zero-sized, so it is not 
clear whether we can trust the perturbative string theory. As we shall see shortly in the 
matrix model, spacetime breaks down and we need to employ a full non-abelian description. 
At later times, strings become strongly coupled, while all dimensions expand, we shall see 
that a rather simple theory emerges, that is, only abelian degrees of freedom survive. 

For case 2, strings are strongly coupled at big bang, transverse dimensions start 
with a zero size, while transverse dimensions x"" have infinite size and contract as time 
evolves. This picture in valid both in terms of the string metric as well as the lOD Einstein 
metric. We will study in more detail the spacetime properties in the matrix model later. 

The string frame metric (^2]) is in general non-flat. f3 never vanishes, so to get a 
flat metric, (i = 1 is necessary, and in this case one chooses (3 = —2a, this is the special 
case considered in [|1|. More generally, one obtains a world-sheet theory with an action 
explicitly depending on time. One can attempt to quantize the string in the light-cone 
gauge, to do so, introduce the new light-cone coordinate dy^ = exp((2a -|- (3)x~^)dx~^ , the 
first term in the string frame metric (^2]) becomes —2dy^dx~ . The lignt-cone momentum 
p_ = —p^ conjugate to x~ is conserved. In the light-cone quantization, we can set the 
light-cone gauge = the bosonic part of the world-sheet action reads 

S = [ drda\e^^''+^^''^d"X'dc,X' + e^^^^a^X^^aX"], (14) 

Ana J 

where the period of a is 27ra'p"'", and exp((2a -|-/3)x+) = (2a-f/3)r. As x~^ — > — oo, r — 0. 
Since r starts at a finite point, it is more useful to use the old light-cone coordinate x~^ = t, 
the action becomes 

Aiva' J 

The time-dependent coefficients may be interpreted as the effective tensions. For transverse 
coordinates X*, the effective tension is T, = -rr^e*-^""'"^^-'*, for coordinates X", the effective 
tension is Ta = 2;^e*^^""'"*^)*. In case 1, both effective tensions get to zero at big bang, 
while the string coupling also gets to zero. Since the string spectrum becomes very dense, 
the usual free string picture does not apply. In later times, effective string tensions get 
large, and the string coupling constant also becomes strong, we will have to deal with a 
strongly coupled massless sector. In case 2, 4a + 2/3 > except for d = 1, our above 
analysis still applies to T^. 2a + 4(3 = 6a(l — 2/^/d), for d > 4, Ta becomes small in 
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earlier time, the same as case 1. For d < 4, Ta becomes large in earlier time, thus these 
transeverse dimensions are effectively frozen, (i = 4 is a special case, is independent of 
time. 

The vertex operator of a string state in the background (|1^) in general is quite com- 
plicated. For instance, consider a massless scalar satisfying equation of motion 

9/. {e-^'^V^g'^'d,^) = 0. (16) 

As function of x"*", the scaling of components is different from the scaling of components 
and g"^^ , so there is no simple plane wave solution in general, which implies that the 
vertex operator of this massless scalar field is not simple. However, we can consider special 
cases when $ is independent of a;", in this vertex operator V = exp(z/c+x"'" + 

ik-X~ + ikix'^) must satisfy the on-shell condition: 

k_{2k+~i-f) - = 0, 7= (9-d)a + d/3, (17) 

which require a imaginary part of k^: Qk- =7/2, thus in the vertex operator, there is an 
exponential factor 

6"^^=^^. (18) 

Since 7 = (9±3v^)a > 0, this exponential factor always blows up at = —00 for d < 9. 
The effective string coupling constant for these states is geff = (/s exp(— ^7^+). Dimen- 
sions for these states are effectively compactified, thus when we discuss interactions 
among these states, the space-time dimensionality is 11 — d, and there is an effective string 
tension in action ([T3|), or a'e// ~ a'exp(— (4a -|- 2(3)x^). The effective Newton constant is 
Geff = gif fOt''^^~'^^^'^ , as a function of x"*", it scales 

Now, 7 is positive and (9 — d){2a -f (3) is non-negative if (i < 9, the effective Newton 
constant blows up at x"*" = —00, the perturbative sstring picture is not valid at least for 
these states. 

Next, we study the matrix model. It is not necessary to compactify any dimension of 
(|]). In the usual flat background, the Matrix Theory action reads 



s = j dm(^^{DtX'f + -^[x\x^Y + ie'^Dte-Re'^-i,[x' 



6] , (20) 



where R is the longitudinal cut-off, or it may be viewed as the radius of x~ in the DLCQ 
M theory For simplicity, we set the M theory Planck length Ip = 1. A derivation of this 
matrix action is given in (see also [|10[ ). With our metric (|^), it is straightforward to 



write the corresponding matrix action. To do this, we need to use the light-cone coordinate 
in (p!0|), and identify with r, the world-line time of DO-branes. In the matrix model, 
we apply the M theory metric (0) directly. Since the action is rather lengthy, we separate 
the action into the bosonic part and the fermionic part. The bosonic part reads 

Sb= [ dTTri^e^^'^^iDrX'f + ^e''(^'=^{D^X''f + ^e^''^^[X\X^]^ 

J 2it 2R 4 ^2]^^ 

4 2 

Note that appearing in the above action is the old light-cone coordinate. The fermionic 
part reads 



Sf^ J dT{i6^Dre - i?e""''6'%[X\ 9] - i?e^"^6''^7a[X", 6]}. 



(22) 



It is rather awkward to use r as time, since it has a finite beginning. Let us switch 
back to the coordinate and on the world-line identify t — x~^ , thus dr = exp{2at)dt, 
we have 

Sb= [ dm{^{DtX^f + J-e2(^-«)*(AX")2 + ^e'^^'iX^X^]^ 

J 2R 2R 4 ^23) 

4 2 

and 

Sf = J dt{ie^Dte - Re^"^9^-f,[X\ 9] - i2e(2«+^)*^^7„[X", 9]}. (24) 

Before study the simplest properties of this matrix model, let us show that we can 
recover the matrix string action of In this case, there is only one X"-, call it X^^. 
The compactification scheme is given in [rT|. Up a dimensionful parameter, we replace 



the trace in the matrix action by J daTv, the commutator R[X^^ , X"^] is replaced by the 
covariant derivative iD^jX^ . Finally, use (3 = —2a, we find 

Sb = j dtdaTv{^{D^X^f + ^glFl -f ^g:^[X\X^]''}, (25) 



and 



Sf = j dtdaTr[9^a''Dc,9 - Rg-^9^-fi[X\9], (26) 



where Qs is the time-dependent string couphng Qs = exp(— 3at). Our action is identical 
to that in [|l| up to an identification of a dimensionful parameter Ig. This matrix string 
theory can be regarded as a 2D Yang-MiUs theory with a time-dependent couphng or a 
2D Yang-MiUs theory with a constant couphng in the world-sheet metric 

ds^ =gj^{-dt^ + da^). (27) 



To study the properties of matrix model defined by ( |23|) and (|2^, we consider cases 
1 and 2 separately. 

• Case 1, /3 = (1 + -^)a. 

The kinetic term of is always simple, as is the kinetic term of 9. The coefficients of 
the remaining terms all vanish in the limit t ^ — oo, so there is no constraint arising from 
these terms, this implies that all matrices are fully non-abelian. On the other hand, as t ^ 
oo, these coefficients blow up, thus, all matrices must commute with one another, the only 
surviving degrees of freedom are diagonal elements. Moreover, X"' must be independent of 
time in this limit, so X"" become frozen abelian moduli. Recall that in the string picture, 
we found previously that the effective tension Ta becomes infinitely heavy, even heavier 
than Ti for large t, this is related to the fact that X" become moduli in the matrix model, 
while X* still have dynamics. If we compactify some of the transverse dimensions, the 
story becomes slightly more involved. For instance, compactifying X^^ on a circle, all 
other matrices become function on a circle a. In the limit t oo, although they have to 
be diagonal, they are not always periodic functions of a, the eigen-values can get permuted 
after circling along cr, these twisted sectors describe strings of various lengths. 

• Case 2, /3 = (1 - -^)a. 

6a and 4a -1-2/3 are always positive, at big bang, there is no constraint on commutators 
[X^,X^], [X^X"], 6'] and [X'^,6]. On the other hand, as t ^ oo, these commutators 
are forced to vanish. 2(/3 — a) is always negative, so at bing bang, X" must be independent 
of time, they become non-abelian moduli in the model. 2a — 4f3 = 6a(l — 2/\/d), d = 4 
becomes the critical dimension. For a larger d, there is no constraint on the commutators 
[X", X^] at big bang, and they have to vanish as t — oo, thus, for d > 4, all matrices 
commute in this limit and only abelian degrees of freedom survive. For d < 4, [X", X^] 
have to vanish at big bang, together with the fact that X" are independent of time, these 
matrices become abelian moduli at big bang, this fact is also refiected in our previous 
analysis in the string picture, where we found become infnitely heavy at big bang. 
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When d = 4, X"' remain nonabelian all the time. In the string picture, recall that when 
d = 4, the effective tension Ta is constant. 

It is also of interest to study compactification. Just as in the flat background, we do not 
know how to write down matrix model action if we compactify more than 5 dimensions. 
We may choose to compactify all X", or some of X", or some of X"" and some of X*. 
For simplicity, let us compactify all X" on a torus T'^. Replacing X'] by iDaX\ 

X'^] by —Fab etc, the bosonic part of the matrix model action reads 

Sb = I d'^+'a{^{Dar - ^e(^"+^^)*(I^„X^)^ + ^.e'^^-'^^^F^^ 



(28) 

The analysis of this action is similar to our previous analysis of the action without com- 



pactification, a statement there can be simply translated to one for the action (28). For 
instance, demanding [X", X^] = is translated to Fab = 0, that is, the spatial connection 
must be flat. 

It is interesting that the previous noticed "critical dimension" d = 4 corresponds to the 
situation that the action (^8[) is no longer complete. Compactiflcation on a torus T^, the 
complete matrix model is to be given by the world volume theory of coincident M5-branes 
|T^. Of course, that theory will also be time-dependent, and hopefully one can flgure out 
more details of the time dependence by studying held theory behavior of (|28|). Likewise, 
action (^8|) is again incomplete for compactification on T^, in this case the complete theory 
is given by the little string theory [|13[ . 

Unlike the case when d = 1, the Yang- Mills theory ( |28| ) can not be simply interpreted 
as a theory with a constant coupling on a world volume of a nontrivial metric. To interpret 
action (P^ ) as Yang-Mills theory, we need to introduce both a nontrivial metric on the world 
volume as well as a time-dependent Yang-Mills coupling. Since a factor gy'ti appears in 
the coefficient of and a factor Qym appears in the coefficient of [X\ X^]"^, we find for 
d> 1 

dsl^y = -e^(4a+2/3)t^^2 ^ ^^(4a+2/3)t^^^a)2^ 
^2^ = e-^(4«+2/3)t+2(a-/3)t^ ^^^^ 

where the subscript WV stands for world volume, the world volume metric components 
Qaa must be all same, there are only 3 independent functions in the geometry and the 
Yang-Mills coupling, however, we need to match 5 different kinds of coefficients in (|28|), 
so it is nontrivial that there is a solution as presented in (|29|). This is quite important 
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especially for d > 3, since we need to complete the action ( pS]) by introducing more degrees 
of freedom, this is possible only when (^8[) admits Yang- Mills theory interpretation. One 
can check that the identification (^9]) also works for the fermionic part of the matrix action. 

Finally, we discuss generalizations of background (^. There are two directions to 
generalize (0). The first is to consider 

ds^ = -2e2"^^dx+da:- + J]e2^'^^(da;^)2. (30) 

i 

The Einstein equations are solved provided 

+ J]A(2a-A) =0. (31) 

i 

This background again preserves 16 supersymmetries of the form e = exp(^ax^)r], with a 
constant ry satisfying 7+77 = 0. 

The second direction to generalize (|^) is to consider metric of the form 

ds^ = e''f^'=^\-2dx+dx- + {dx')'') + e2'?("^)(rfa;")2. (32) 

The non- vanishing curvature 2 forms are 

(33) 

Ra+ = i2f'g'-g''-g")e-'fe'^Ae+. 
The Einstein equations are solved if 

(9 - d){f" - n + d{2f'g' - g'' - g") = 0. (34) 

Once again, there are 16 unbroken supersymmetries, parameterized by e = exp(i/)r7, 
7+77 = 0. 

Needless to say that equation (|3^) admits infinite many solutions. Two special solu- 
tions deserve attention. One solution is when g = 0, in this case /'^ — /" = 0. All the 
curvature 2 forms vanish thus it may appear that we obtain a fiat space solution. This is 
almost case except for a simple twist. The solution of f'^ — f" = is exp(2/) = l/(x+)^, 
by changing coordinates dy^ = dx~^ / {x~^)^ , y~ = x~ we obtain a metric 



ds'^ = -2dy+dy- + {y+f{dx'f + {dx^f. 
9 



(35) 



If one of x'^ is periodic, we find a nontrivial background similar to the orbifold discussed 
in 0,0. Of course if all are noncompact, (^5]) is equivalent to the Minkowski space. 
Similarly, one may set / = and find exp(25f) = (x"'")^, a nontrivial orbifold results 
provided one of a;" is compactified. One may attempt to construct a matrix model for 



the metric ( |35|) too. However, time t = seems to have to terminate at r = 0. One 
may try to avoid this problem by going back to the original light-cone coordinate and 
take t = in the matrix model, although this helps to eliminate the coefficient in the 
kinetic term of X* but it also introduces a coefficient in the kinetic term of X", again t 
must be terminated at t = 0. Also the supersymmetry parameter e = \J x'^r\ does not exist 
beyond x'^ = since e has to be real. This problem may cause disease in the definition of 
the matrix model, and may be a reflection of the gravitational instability caused by a test 
particle. 
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